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Abstract. The thermal, spin and magnetic evolution of neutron stars in the old low 
mass binaries is first explored. Recycled to very short periods via accretion torques, the 
neutron stars lose their magnetism progressively. If accretion proceeds undisturbed for 
100 Myrs these stars can rotate close to break up with periods far below the minimum 
observed of 1.558 ms. We investigate their histories using population synthesis models 
to show that a tail should exist in the period distribution below 1.558 ms. The search 
of these ultrafastly spinning neutron stars as pulsars can help discriminating among 
the various equations of state for nuclear matter, and can shed light into the physics 
of binary evolution. 

The evolution of isolated neutron stars in the Galaxy is explored beyond the pulsar 
phase. Moving through the tenuous interstellar medium, these old solitary neutron stars 
lose their rotational energy. Whether also their magnetism fades is still a mystery. A 
population synthesis model has revealed that only a tiny fraction of them is able to 
accrete from the interstellar medium, shining in the X-rays. There is the hope that 
these solitary stars will eventually appear as faint sources in the Chandra sky survey. 
This might give insight on the long term evolution of the magnetic field in isolated 
objects. 



1 Introduction 

The amount of rotation and the strength of the magnetic field determine many 
of the neutron star's observational properties. Over the neutron star lifetime, the 
spin and the field change and the study of their evolution provides important 
clues into the physics of the stellar interior. 

Of the billion neutron stars in the Galaxy, we shall be mainly concerned 
with the evolution of two distinct populations: The millisecond pulsars and the 
isolated neutron stars both aging, the first in low mass binaries, the second as 
field stars moving in the interstellar medium of our Milky Way. We will show 
that the interaction with their surroundings may profoundly alter their spin 
and magnetic field. The extent of these changes and the modes vary in the two 
scenarios. It is in exploring this diversity that we wish to infer the nature of 
the equation of state and to provide a unified view of field decay. In particular, 
the existence of "unconventional" sources, such as sub-millisecond pulsars, stars 
rotating close to their break up limit, and as solitary neutron stars accreting 
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the interstellar medium, is a crucial test for our studies. Their discovery is an 
observational and theoretical challenge. 

In Section 2 we trace the evolution of a neutron star in the Ejector, Propeller 
and Accretion phases, described briefly using simple background arguments. 
Section 3 surveys the physical models for the evolution of the magnetic field both 
in isolated and accreting systems. In Section 4 we explore possible individual 
pathways that may lead to the formation, in low mass binaries, of neutron stars 
spinning very close to their mass shedding threshold, a limit sensitive to the 
equation of state for nuclear matter. Pathways of isolated neutron stars follow. 
The first four sections set the general background used to construct, in Sections 
5 and 6, statistical models aimed at determining the presence and abundance of 
sub-millisecond pulsars in the Galaxy, and of solitary neutron stars shining in 
the X-rays. Specifically, in Section 5 we explore, within the recycling scenario, 
the star's spin and magnetic evolution using physical models and the role played 
by disc instabilities in affecting the latest phases of binary evolution toward the 
Ejector state of sub-millisecond pulsars. In Section 6 we carry out the first stellar 
census. We then establish how elusive neutron stars can be as accreting sources 
from the interstellar medium due to their large velocities and to magnetic field 
decay. 

2 The Ejector, Propeller and Accretion Phases 

Over the stellar lifetime, magnetic and hydrodynamic torques acting on the 
neutron star (NS hereafter) induce secular changes in its spin rate. Four physical 
parameters determine the extent of the torques: the magnetic field strength B, 
the rotational period P, the density of the surrounding (interstellar) medium 
n, and the rate of mass inflow M toward the NS. According to the magnitude 
of these quantities, a NS experiences three different evolutionary paths: Ejector 
(£), Propeller (V), and Accretion (.4) In phase £ the NS braking results 
from the loss of magneto-dipole radiation as in an ordinary pulsar. The implied 
spin-down rate is 



where / is the NS moment of inertia and fi = B s R^/2 the NS magnetic mo- 
ment, function of the stellar radius i? s and polar magnetic field B s . The torque 
decreases with increasing P, and at the current period P (much longer than the 
initial period) the time spent in phase £ is 



dP 
dt 
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Resulting from the emission of electromagnetic waves and charged particles, the 
rotational energy loss of equation (1) can proceed well beyond the active radio- 
pulsar phase and the magneto-dipolar outflow creates a hollow cavern nesting 
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the NS. Phase £ remains active as long as the characteristic radius r st of this 
cavern (the stopping radius || B) is larger than both the light cylinder radius 



cP 



and the gravitational radius 



2GM 

r G = 2 , 2 4 

where c s is the sound speed of the interstellar medium (ISM) and v Te \ the NS 
velocity relative to the medium. Phase £ ends when matter leaks through r st , 
and this occurs when the period exceeds 

p - - ^" ,/v/ V' /4 < 5 ' 

where v = (c 2 + v 2 ^) 1 / 2 and g is the mass density of the ISM 0). At Ps^-p 
the outflowing pulsar momentum flux is unable to balance the ram pressure 
of the gravitationally captured matter at re 0. Matter initiates its infall and 
approaches the magneto spheric radius (or Alfven radius) 

V 2GMM 2 J 

Here the magnetic pressure of the dipolar stellar field (oc r~ 6 ) balances the ram 
pressure of the infalling material (oc r~ 5 / 2 ), accreting at a rate M. At ta the 
steeply rising magnetic field would thread the flow, enforcing it to corotation. 
However, not always can matter pass this edge: penetration is prevented when- 
ever the rotational speed Q of the uniformly rotating magnetosphere exceeds the 
local Keplerian velocity ujk (= [GAf/r^] 1 / 2 ) at r&. This condition translates into 
a comparison between the corotation radius 



(GM\ 



1/3 

(7) 



and r\. When r cor lies inside ta the magnetosphere centrifugally lifts the plasma 
above its escape velocity, inhibiting its further infall toward the NS surface: this 
is the propeller phase V . The magnetically driven torques lead to a secular spin- 
down of the NS occurring at a rate 

j t (W) = CMr c 2 01> K(r A ) _ Q]j ( 8 ) 

where M is the accretion rate at ta, and £ a numerical factor dependent on 
the accretion pattern, ranging from £ ~ 1 for disc accretion to £ ~ 10~ 2 for 
spherical accretion. The effectiveness of the propeller in spinning down the NS 
and in ejecting matter far out from the magnetosphere is still largely model 
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dependent Q. Recently, evidence that this mechanism is at work in binaries of 
high Q and low mass j(| has come from X-ray observations. 

Phase V terminates only when the corotation radius r cor increases above 
ta; thereon the NS is able to accept matter directly onto its surface: this is the 
accretion phase A. The transition V — > A occurs at a critical NS spin rate P-p^A 
obtained by equating r A and r cor : 



V^A 



2 * .(-SJT* *!L. (9 ) 
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In isolated NS, the accretion flow is almost spherically symmetric and the 
NS period P, after the onset of A, may vary erratically due to the positive and 
negative random torques resulting from the turbulent ISM @ § 

In a disc-like geometry (in binaries) the spin evolution during accretion is 
guided by the advective, magnetic and viscous torques present in the disc. De- 
pending on their relative magnitude, these torques can either spin the NS up or 
down. The rate of change of the NS angular momentum is obtained integrating 
the different contributions over the whole disc, yielding 

^-(W) ~ Mr 2 A cu K (r A ) + / B 2 — [u K (R) - Sl]dR , (10) 
at J rA r\ 

where r\ is the plasma electrical diffusivity, h the disc scale thickness and R the 



radial cylindrical coordinate (in eq. 10 we neglect the local viscous contribu- 



tion as it vanishes at r A ). The first term in the rhs of equation ( [LOP describes 
the advection of angular momentum by accretion. Magnetic torques are instead 
non-local, resulting from an integral that extends over the whole disc: different 
regions can give either positive (spin- up) or negative contributions (spin-down). 
When SI lok^a), the integrand in ([l0]) is negative yielding a spin-down mag- 
netic torque. On the other hand, if the star rotates very slowly (SI <C ojk{ta)) 
the overall magnetic torque is positive, leading to a secular spin-up: when the 
total torque is positive the spin-up process is termed recycling. At a critical (in- 
termediate) value of SI, termed equilibrium angular velocity u> oq , the magnetic 
torque is negative, and its value exactly offsets the (positive) advective torque. 
In this regime the NS accretes matter from the disc without exchanging angular 
momentum at all. The precise value of co eq is rather uncertain, but it ought to 
lie between 0.71 — 0.95 times the critical frequency cj cr = 2tt / P-p^^, in the mid 
of the propeller and accretor regimes (see |9J and references therein) . Since ui CT 
and uj cq have near values, they are often confused in the literature. There is a 
narrow range for SI between w oq and w cr within which the NS accretes mass but 
spins down. 

Once the star has reached its equilibrium period P oq = 2n/ij cq , further 
changes in the spin period P are possible only if the magnetic field and/or 



1 Old accreting isolated NSs are thus expected to show strong fluctuations of P over 
a time scale comparable to the crossing time rfl uct ~ min(nsM, rG)/v rc i, where hsm 
is the spatial scale of the ISM inhomogeneities. 
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M vary. Generally, field decay causes the star to slide along the so called spin-up 
line where the period P takes the value P cq (n, M), at the current values of /i and 
M. A large increase in the mass transfer rate can induce further spin-up since 
the magnetosphere at ta is squeezed in a region of higher Keplerian velocity. 
A decrease in M has the opposite effect and if there is a decline in the mass 
transfer rate (due to disc instabilities and/or perturbations in the atmosphere of 
the donor star), the NS can transit from A — > V — > £. The critical period for the 
last transition does not coincide with equation (0), but occurs when the Alfven 
radius ta (function of the mass transfer rate M) exceeds the light cylinder radius 
r\ c (eqs. |^ and i.e., when 

c \2GMM 2 ) V ' 

(note that this transition is not symmetric relative to £ — > V). Over the stellar 
lifetime, NSs can trace loops moving through the various phases (V <-> A or/and 
V «-» £). It may happen that ta increases above tq when M is very low (as in 
the ISM): in this case the star is in the so called Georotator state (Q). 

The major drivers of the NS evolution are magnetic field decay (though 
amplification is an interesting possibility) and variations in the mass transfer 
rate: their study is thus the subject of §3 and §4. 



3 Magnetic field evolution 

In the /it versus P diagram of pulsars (PSRs) we clearly find NSs endowed with 
intense magnetic fields (the canonical pulsars) and NSs with a much lower mag- 
netic moment, the "millisecond" pulsars (MSPs). The first are isolated objects 
and largely outnumber the millisecond pulsars often living in binaries with de- 
generate companion stars. The origin of this observational dichotomy seems un- 
derstood but the physical mechanisms driving the field evolution remain still 
uncertain. We review here current ideas for the evolution of the B field. They 
are applied later when investigating (i) the existence of sub-millisecond pulsars, 
and (ii) the nature of the six isolated NSs discovered by Rosat. 



3.1 Historical and observational outline 



Soon after the discovery of PSRs, Ostriker & Gunn 10 proposed that their sur- 
face magnetic field should not remain constant but decay. Under this assumption 
they explained the absence of PSRs with periods much longer than one second 
(the decay of /i would bring the objects below the death line, turning off the 
active radio emission |ll[] ) and the claimed dimming of the radio luminosity L# 
proportional to /z 2 j[2) . Their pioneering statistical study suggested that the NSs 
at birth have very high magnetic moments fj, ~ 10 29 -j- 10 31 G cm 3 , whose values 
decay exponentially due to ohmic dissipation in the NS interior, on a typical 
time-scale lO] 
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T„ = 4ffi? s 2 /7TC 2 , (12) 

where a is the conductivity, taken as uniform. Ostriker & Gunn evaluated r M 
using the electrical conductivity for the crystallized crust (~ 0.6 x 10 23 sec -1 
[fl4|), thereby assuming implicitly that the magnetic field resides in the outer 
layers (p < 10 14 g cm -3 ) of the star, and obtained r M ~ 4 x 10 6 yr. But in the 
same year, Baym, Pethick & Pines Q argued that the magnetic field pervades 
the entire star and so one had to adopt the much higher a of the core, thus 
predicting a very long r M exceeding the Hubble time, a result confirmed ]l6j in 
1971 when different conductivities in the NS interior were taken into account. 
When a larger number of pulsars became available for statistics, many authors 
compared the observations with population synthesis Monte Carlo models. At 
first data seemed to support evidence for field decay, but subsequent inves- 
tigations |Q Jl9| |2(J Q have indicated that the observed PSR population is 
compatible with no decay or with decay on time-scales r M > 10 8 yr longer than 
the pulsar phase, beyond which the NSs become practically invisible. There is 
a way to probe indirectly the possible decay of the field by searching for those 
old NSs moving with low speed (< 40 km s _1 ) that can gravitationally capture 
the interstellar medium and shine as dim accreting sources. As discussed in §6, 
their presence in the Galaxy depends on the evolution of the field. 

As far as the old MSPs are concerned, these sources possess very low values 
of fi ~ 7 x 10 25 — 10 27 G cm 3 ||^] when they appear as pulsars. The observations 
are again consistent with no field decay over their radio active phase, a result 
inferred from the old age of the dwarf companion stars, that implies a at least 
comparable to their true galactic age Q Q [||| ||(| . This does not preclude from 
the possibility that field decay has occurred during their previous evolution: this 
is discussed and described in §3.4. 

In the next sections we survey the models for the decay of the magnetic field 
both for isolated NSs and for the NSs in binaries. 

3.2 Field evolution in isolated objects: spontaneous decay? 

Models for magnetic field decay in isolated NSs flourished over the last decade. 
A first improvement was to abandon the hypothesis of a uniform conductivity in 
the stellar crust and core; a increases when moving from the outer liquid layers 
to the deeper crystallized crust J27|| . Sang & Chanmugam [^8| noticed that the 
decay of a field residing only in the crust is not strictly exponential because of its 
inward diffusion toward regions of higher a. Moreover, the decay rate depends 
on the (highly uncertain) depth penetrated by the initial field: the decay is more 
rapid if, at birth, the field resides in the outer lower density regions of the crust 
both because a is lower and because t m depends on the length scale of field 
gradients (less than a tenth of the stellar radius R if fi is confined in the crust). 

The inclusion, in the calculation, of the cooling history of the NS was a second 
major improvement |29]| ]3C[| | |3l) . It confirmed the non exponential behavior 
of the decay, showing in addition that a slower cooling would accelerate the 
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decay (as a warmer star has a lower conductivity). A first short phase ( ;$ 10 6 
yr) of comparatively rapid decay was found, during which \x drops from values 
~ 6 -j- 30 x 10 30 G cm 3 , typical of young pulsars in supernova remnants, to more 
canonical values ~ 1-^3 x 10 30 G cm 3 typical of normal PSRs. (During this phase 
the electron-phonon interaction dominates a) . A second phase of no decay was 
found lasting 10 -j- 300 Myr, followed by a power-law decay (dominated by the 
interaction of the electron upon the impurities) and eventually an exponential 
decay. 

If the field penetrates the whole star, the core ohmic decay time-scale would 
be too long to allow decay p2| . Thus, an alternative scenario was proposed on 
the hypothesis that NSs have superfluid and superconducting cores. The angular 
momentum of the core is believed to be carried by |33| A vort ex — 2 x 10 16 /P sec 
neutron vortex lines, parallel to the spin axis; the magnetic flux should be con- 
fined in A^uxoid quantized proton flux tubes, where Afl uxo ;d — 10 32 Bc,u, with 
Bc,i3 the average field strength in the core, expressed in units of 10 13 G. In 
1990, Srinivasan et al. Q suggested that the inter-pinning Q between these 
quantized entities causes the fluxoids to be carried and deposited in the crust 
as the NS spins down. The magnetic flux penetrated into the crust will then 
decay due to ohmic dissipation and the surface magnetic field of the NS keeps 
decreasing until it relaxes at the value of the residual field in the core. In this 
framework the spin history of the star drives and controls the changes of the 
magnetic field. Assuming the same radial velocity for both the fluxoids and the 
vortices this model implies a time-evolution 

n(t) = no ( i + t/ Tsdd y 1/4 (13) 

occurring on a spin down time 

r sdd ~ 8 x 10 6 P 2 /45/V30 yr, (14) 

where index denotes values at the onset of phase £, and I45 the NS moment of 
inertia in units of 10 45 g cm 2 . Equation ( [l3"| ) models the simplest version of the 
spin down driven (sdd) decay of /z in isolated NSs. It can be refined, accounting 
for the non instantaneous relaxation of the surface magnetic field to the value 
in the core Q. 

Disappointedly, until now it has not been possible to perform a reliable statis- 
tical test of all the models. They predict clearly distinct evolutionary pathways 
only during the early stage of the pulsar life, namely during the first 10 6 yr after 
birth. Therefore a large sample of young NSs is requested for a comparison. Up 
to now, the pulsar catalog lists only an handful of such objects, as the observed 
PSR population is dominated by elder sources, with characteristic time of few 
10 7 yr. In such situation, there is no statistical clue for rejecting the hypothesis 
of a non decaying field. Some of the surveys in progress, as the one running at 
Parkes p7| , dedicated to the search of young PSRs, could ultimately solve this 
problem. 
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3.3 Magnetic field evolution in binaries: accretion driven scenario 

None of the models for spontaneous decay, proposed so far for the isolated ob- 
jects, can explain magnetic moments of ~ 10 26 G cm 3 , which are typical of the 
millisecond PSRs. In order to attain these values of /i over a Hubble time, the 
crustal models would require currents located only in a narrow layer below a 
density p < 6 x 10 10 gem -3 , a very unlikely possibility. On the other hand, the 
"sdd" models would predict only a modest decrease of /x, at most of two orders 
of magnitude with respect to the initial values. However, one common feature 
groups the low-/i objects, i.e., that all the PSRs with low magnetic moment spent 
a part of their life in a multiple stellar system. Guided by this fact, many authors 
have tried to relate the low values of /j, with the interaction between the NS and 
its companion star(s). Two different scenarios have been proposed so far: (i) The 
accretion driven scenario developed under the hypothesis of a crustal magnetic 
field; (ii) The spin driven scenario lying on the hypothesis of a core magnetic 
field. 

In this section we consider scenario (i) taking the hypothesis that the cur- 
rents generating the field are located deep in the crust []. A first observational 
suggestion about the effects of the accretion on a (crustal) magnetic field was 
given in 1986 by Taam & van den Heuvel [Q: they noticed an approximate 
dependence of the surface magnetic field on the amount of mass accreted onto 
the neutron star. Shibazaki et al. J4l| ] later presented an empirical formula for 
the decay in presence of accretion: 

_ IM) (15") 

1 + AMacc/m* 

where po represents the magnetic moment at the onset of accretion, AM acc is the 
amount of accreted mass and a parameter to be fitted with the observations. 
These authors claimed that m* ~ 10~ 4 M Q could reproduce the Taam & van den 
Heuvel's correlation. Zhang, Wu & Yang (4^] [Q gave physical foundation to ( p"5| ) 
assuming that the compressed accreted matter has ferromagnetic permeability. 
However, using a larger database Wijers Q] showed that these models are not 
fully consistent with the available data both on X-ray binaries and recycled 
pulsars. Romani |^5| first introduced the accretion rate M as a second parameter 
in driving the magnetic moment decay in addition to AM acc . He pointed out that 
the accretion produces two major effects: (i) heating of the crust (depending on 
M), which determines a reduction of the conductivity (and in turn the hastening 
of ohmic decay) and (ii) advection of the field screened by the diamagnetic 
accreted material. Moreover, the advection of the field lines stops when /i <, 10 27 
G cm 3 , resulting in an asymptotic value for the magnetic moment (in agreement 
with the observations). 



The onset of a thermo-magnetic instability, which transforms heat into magnetic 
energy at the moment of NS formation, is an effective mechanism to produce strong 
fields in the crust of a NS |^] [Q. Although this instability is not yet completely 
explored for poloidal fields, it is a plausible mechanism for the origin of a crust field 
which does not depend on special assumptions. 
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The existence of a bottom field (with a predicted scaling fi oc M 1 / 2 ) is 
also possible if the currents sustaining the surface magnetic field of the NS 
are neutralized by the currents developed in the diamagnetic blobs of accreted 
matter In this case the decay of p typically ends when the accretion disc 
skims the NS surface p7[ . 

The first fully consistent theoretical calculations of the accretion-induced- 
decay of the crustal magnetic field of a NS was performed by Urpin, Geppert & 
Konenkov |h| |Q |HJ . The basic physical mechanism is the diffusion (ohmic 
decay of currents) and advection of the magnetic field sustained by currents 
circulating in the non superconducting crust. The magnetic field evolution is 
calculated according to the induction equation: 

a td 2 / I \ 

-VA -VAB UVA(bAB) (16) 



dt 4tt 



a 



where B is the magnetic induction and v is the velocity of the moving fluid (= 
in a non accreting NS, \v\ = M/4irr 2 p for a radial approximation of the accreted 
fluid flow). In their calculation, the superconducting NS core is assumed to expel 
the magnetic field (in the following we will refer to this boundary condition for 
the field at the crust-core interface as BCI); the induction equation is solved for 
a dipolar field and a is given as in j5S| . 

Accretion affects the ohmic decay in two ways: (i) heating the crust, so re- 
ducing a; (ii) transporting matter and currents toward the core-crust boundary. 
The relative importance of these two effects depends on other physical parame- 
ters, such as the initial depth penetrated by the currents, the impurity content 
(Q), the accretion rate, the duration of the accretion phase and the equation of 
state of the nuclear matter. In coupling the magnetic history of the star to its 
spin evolution (through £ to A), Urpin, Geppert & Konenkov Q found that 
neutron stars born with standard magnetic moments and spin periods can evolve 
to low-field rapidly rotating objects, as illustrated in Figure [l| 

The behavior of currents and the response of nuclear matter to an interior 
field, at the interface between the crust and the core of a NS (typical density 
~ 2 x 10 14 g cm -3 ), is still an open issue for the theorists. In view of these 
uncertainties, Konar & Bhattacharya J54| chose a different boundary condition 
at the crust-core interface for solving the equation (0). They noticed that the 
deposition of accreted matter on top of the crust can imply the assimilation of 
original current-carrying crustal layers into the superconducting core of the NS. 
In particular they assumed J54j] the newly formed superconducting material to 
retain the magnetic flux in the form of Abrikosov fluxoids Jig] rather than to 
expel it through the Meissner effect (in the following this boundary condition 
will be labeled as BCII). Within this model, accretion produces a third effect 
on the evolution of /x: (iii) the assimilation of material into the core, where the 
conductivity is huge, freezes the decay. 

As illustrated in Figure ||, also this hypothesis leads to a reduction of the sur- 
face magnetic field of more than 4 orders of magnitude in 10 8 yr, explaining the 
low fi and its possible freezing in the millisecond pulsar population. Likewise for 
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Fig. 1. Magnetic, rotational and thermal evolution vs time of a 1.4 Mq NS with /i = 
10 30 G cm 3 at birth, as calculated in under BCI. The adopted equation of state 
is PS (a stiff one). The initial depth of the currents is 10 13 g cm -3 and the impurity 
parameter Q — 0.01. Phase I and II refers to £ and V respectively. The value of M 
from the stellar wind (referred to as phase III of A) is 2/3 x 10 _8 Me, and the lifetime of 
the companion star on the main sequence is 5 x 10 9 yr. During the disc accretion phase 
(referred to as phase IV), M is 2/3 x 10 _3 M E (solid line) or 2/3 x 10~ 2 M E (dotted 
line). Phase IV lasts 10 s yr. The inserted panel enlarges the spin evolution during the 
very short phase at the boundary between phase III and IV, that is when disc accretion 
sets in. 



BCI, it emerges that the decay of \i depends not only on the total accreted mass, 
but on the accretion rate itself. However, in this case the higher the accretion 
rate is, the stronger is the pull of crustal material into the core resulting in an 
earlier freezing of the decay. In a subsequent paper |^5| , it was claimed that such 
a positive correlation between M and the final value of /i is supported by some 
observational evidence |5(| ]57| . 

Other refinements of the accretion-induced decay scenario include the effects 
of a non spherical symmetric accretion , the evolution of multipoles at the NS 
surface |38| , the post-accretion increase of /i due to the re-diffusion of the buried 
field towards the surface ]5j|, and the relativistic corrections to the decaying 
history All together, they produce only slight changes to the aforementioned 
results. 
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Fig. 2. Evolution, under BCI (left) and BCII (right panel), of the surface magnetic 
moment /j,(t) (relative to its initial value) of a 1.4 Mq NS undergoing accretion at five 
different levels. Top (bottom) curve refers to M = 1CT 11 Mq y~ x (M = 10~ 9 M y" 1 ), 
from pM. Equation of state "PS" is adopted. 



3.4 Magnetic field evolution in binaries: spin driven scenario 

Both the simple idea of the flux conservation during the gravitational collapse 
and the action of a dynamo in the convective proto-neutron star [ |62[ lead to 
the existence of a magnetic field penetrating the whole NS. Due to the huge 
conductivity of the core matter, no ohmic decay occurs during a Hubble time 
and so other effects have been invoked to account for the low field NSs. 

Besides magneto-dipole braking, a NS experiences a longer and more signif- 
icant phase of spin-down, phase V (described in §2) which enhances the effect 



of the spin-down driven scenario sketched in §3.2. This is particularly impor- 
tant when the NS lives in a binary system. Miri & Bhattacharya |36| and Miri 
[ |63| explored the case of low mass systems showing that at the end of V and 
A phases the magnetic moment /i has decayed, relative to its initial value /zo, 
by a factor scaling with Po/Pm&x where Pq is the initial rotational period and 
-fmax the longest period attained during the phase where V is active (V results 
from the interaction of the companion wind with the NS). Assuming that the 
star spins-down to P 10 3 sec, a residual magnetic moment is 



^finai * \\ m Mo * fO- 4 - lO- 3 ^ * 10 26 -j- 1CJ 27 G cm 3 (17) 
1000 sec 

compatible with the values observed for MSP population. 

Recently, Konar & Bhattacharya |64| re-examined this scenario in the case 
of accreting NSs, incorporating the microphysics of the crust and the material 
movements due to the accretion. They concluded that the model can reproduce 
the values of /i observed in the low-mass systems only for large values of the 
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unknown impurity parameter Q £ 0.05 (in contrast with accretion driven decay 
models demanding for values of Q 0.01). If the wind accretion phase is short 
or absent, Q should even exceed unity. A more serious objection to this model 
was risen by Konenkov & Geppert |65) who pointed out that the motion of the 
proton flux tubes in the core leads to a distortion of the field structure near the 
crust-core interface and this in turn creates a back-reaction of the crust on the 
fluxoids expulsion: it results that (i) the sdd can be adequate for describing the 
//-evolution only for weak initial magnetic moments (/l*o 10 29 G cm 3 ) and (ii) 
the predicted correlation //e na i oc -P^ax is no more justified. 

Finally, we have to mention another class of models in which the magnetic 
evolution is strictly related with the spin history of the star: in 1991 Ruderman 
[ |66| J(37j proposed that the coupling among P and /i could take place via crustal 
plate tectonics. The rotational torques acting on the neutron star cause the 
crustal plates to migrate, thus dragging the magnetic poles anchored in them. 
As a consequence the effective magnetic dipole moment can strongly vary, either 
in intensity (decreasing or growing) and in direction, suggesting a tendency to 
produce an overabundance of both orthogonal (spin axis _L magnetic axis) and 
nearly aligned rotators (spin axis |j magnetic axis). Current observation of disc 
population MSPs seem compatible with this feature |||. Moreover the crust- 
cracking events could account for the glitches seen in young PSRs ]69[| . A more 
extended description of the variety of the spin driven decay models can be found 
for example in |7(J. 

4 Evolutionary pathways in various environment 

We here describe shortly possible NS pathways driven by the changes in the field 
and in the accretion rate, over phases £, V or A, in low mass binaries and in the 
ISM. 

4.1 Binaries 

The NSs in binaries experience a complex evolution which is tightly coupled to 
the orbital and internal evolution of the companion star jn| . The powerful pulsar 
wind initially sweeps the stellar wind away and the NS spends its lifetime in the 
£ phase. With the weakening of the electro-magnetic pressure with increasing 
period, an extended V phase establishes (lasting 10 s — 10 9 yr) during which the 
NS is spun down further to periods of ~ 10 2 — 10 4 seconds in its interaction 
with the stellar wind of the companion star. When the period P-p^A is attained, 
accretion sets in down to the NS surface and we could possibly trace this phase 
identifying the X-ray emission from the wind fed NS. Here the NS can accrete 
with no exchange of angular momentum with the incoming matter unless the 
magnetic field decreases due to either accretion induced or spin induced decay. 
The NS can slide along the corresponding equilibrium spin-up line. As soon as 
the donor star evolve into a giant state, it can fill its Roche lobe and matter 
overflows from the inner Lagrangian point forming an accretion disc. When a 
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disc establishes (on the viscous time scale) large accretion rates are available and 
from this moment on recycling starts. In the Roche lobe overflow phase (RLO) 
the NS can be re-accelerated to millisecond periods while the field decays down 
to values of 10 8 — 10 9 G. The mass transfer and orbital evolution are complex to 
model and it is now believed that a significant fraction of the mass lost by the 
donor does not accrete onto the NS, but is ejected from the system JT2|. 




d ■ ■ — r r I — i i rmt i l — i 

1CT a 10- s 10" 1 10° 10 1 10 2 10 3 

P[s] 

Fig. 3. The evolutionary tracks of NSs in the plane B — P for different durations of 
the evolution (labeled by the numbers on the right of the curves, in logio yr) during 
phases I(£) + II(P) + III (A) + IV (A). The initial surface magnetic field is 3 x 10 12 
G, Q = 0.03, and the other parameters as in Figure |l] |5l| . 

Once reached the spin-up line at the current accretion rate, the NS loads 
matter acquiring angular momentum at the rate of field decay. Recycling ceases 
when the donor star has evolved beyond the red giant phase and a dwarf remnant 
is left over stellar evolution, or when the binary as become detached, so that the 
donor underfills its Roche lobe j7l). If accretion ends abruptly, the NS may 
avoid phase V, transiting directly to £ and possibly re- appearing as an active 
radio pulsar. In low mass binaries where the donor star is a low main sequence 
star, evolution proceeds in a fairly ordered way: from £ to V (from the stellar 
wind), to A (from wind fed accretion), to A again (from a Keplerian disc; RLO), 
as shown in Figure |[ Eventually the NS transits to £ or V j£ when accretion 
halts. Observational and theoretical considerations hint for a large decay of the 
magnetic moment [i in phase A, particularly during RLO, as outlined also in 
Figures 1 and 3. The type of evolution will be applied in §5 to address the 
problem of the existence of sub- millisecond pulsars. 
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4.2 Isolated neutron stars 



At birth, isolated NSs experience phase £. The subsequent evolution depends 
on the NS velocity relative to the ISM, the value of the ISM density and the 
magnetic field intensity. The large mean kick velocities acquired at birth ((V) ~ 
300 km s" 1 , Q dlUZllB) can re duce the extent of magnetic torques in phase 
V and can impede accretion fully. It is thus possible that phase £ never ends. 
This occurs when the NS speed V is as high as V > 100 /^on 1 / 2 km s -1 : This 
speed is derived estimating the duration of phase £ at Ps-*-p (eqs. 2 and 5) 



10M30V 



(18) 



imposing r £ ^-p ~ 10 yr and v ~ V (in eq. 18, V\o, n, and /130 are the velocity, 
density, and magnetic moment in units of 10 km s -1 , 1 cm -3 , and 10 30 G cm 3 , 
respectively; hereon we will denote the normalizations as subscripts, for sim- 
plicity). Clearly, only the slowest NSs can loop into the accretion phase under 
typical ISM conditions ( where n ^ 1 cm -3 ). Molecular clouds seem a privileged 
site for phase A J7f| , but the probability of crossing these high density regions 
is relatively low and shortlived. 




8 10 



Fig. 4. Phases £ , V , A and Q, in the Vio — M30 plane for isolate d N Ss moving in a ISM 
of density n = 1 cm -3 and a sound speed of 10 km s _1 , from [ 10 1 1 . 



Typically, in the ISM accretion is nearly spherical and a Bondi-Holye type 
flow establish with 

i2 „ . r 2„ t„1 , v 2 \-3/ 2 



(19) 



m = ■>z(;-.\rf>(c.i -r , lvl , 

where here w rc i coincides with V. The critical periods thus are 

P £ ^ v - \0n- 1/A V^ 2 ny Q 2 s and P V ^A ~ 500 TT 3/r V^ 7 fi% 7 s, (20) 

if the magnetic field is constant in time, c s ~ 10 km s -1 , and M = 1.4 M Q . 
As indicated in equations (18-20), the NS evolution is guided by the two key 
parameters, V and /i (at fixed ISM's density). This is quantitatively illustrated 
in Figure where the NS's phases are identified in the V vs B plane. A strong 
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constant magnetic field (fx > 10 30 G cm 3 ) implies large magneto-dipole losses 
that decelerate rapidly the NS to favor its entrance in phase A (even when 
V ~ 300 km s -1 ), despite the long period of the V — ► A transition (note eq. 20). 

What are the consequences of field decay on the evolution? Whether field 
decay enhances or reduces the probability of a transition beyond £ is a subtle 
question that has recently been partly addressed |T6| |77j JFsj Q . Two competing 
effects come into play when B decays: (i) the spin-down rate slows down because 
of the weakening of B, causing the NS to persist longer in state £; (ii) the periods 
Pg^-p and P-p^A instead decrease with B, and this acts in the opposite sense. 



I 



5 



Bottom magnetic field. G 



Fig. 5. Loci, in the r M — £?bo plane, where phase £ lasts longer than 10 yr. The curves 
refer (from top to bottom) to initial fields of 5 x 10 11 (dotted line), 10 12 (dashed line), 
2 x 10 12 G (hatched region), for v IC \ + c a — 40 km s" 1 , and ISM density n = 1 cm' 3 . 



To explore this delicate interplay, Popov & Prokhorov [Q studied a toy 
model where a field is exponentially decaying on a scale r M from its initial value 
/Zij down to a bottom value ^bo- Under this hypothesis, the Ejector time scale 
(with a decaying fi) can be estimated analytically to give 




= 2r bo /r t 



TS.a > T bo 



(21) 



where T bo = T,jln(/io//Xbo) is the time when the bottom field is attained and T£ 
is from equation (2) for a constant field. Figure shows the loci where the time 
spent in phase £ equals the age of the Galaxy, i.e., 10 Gyr, as a function of the 
bottom field B bo and of the decay time-scale r M . The "forbidden" region lies just 
below each curve. As it appears from the Figure, the interval over which the NS 
never leaves stage £ is non negligible: It widens when /xp is decreased, due to 
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the weakness of the magneto-dipole torque. Large as well as very low bottom 
fields do not constrain and permit entrance to phase A. When /ibo is very 
low, there is a turn-off of all magnetospheric effects on the inflowing matter and 
matter accretes promptly. 

Similarly, Colpi et al. |7(| considered a model in which the spin evolution 
causes the core field to migrate to the crust where dissipation processes drive 
ohmic field decay (Sect. 3). In this circumstance, the entrance to phase A is less 
likely if r M ~ 10 8 yr, but possible otherwise. Thus, field decay can hinder the 
stars in phase £ if its decay is somewhat "fine-tuned", while a fast decay would 
drive them into A promptly. The strongest theoretical argument against phase 
A remains however the high kick velocity that the NSs acquire at birth (see the 
review of Lai in this book) . There remain nevertheless open the possibility that 
weak field accreting NSs exist in the Sun's vicinity. 

We now proceed on exploring the existence of "unusal" NSs, the sub-MSPs 
and the X-ray isolated NSs for which these pathways are statistically relevant. 

5 A population synthesis approach to the formation of 

sufe-MSPs 

This section is devoted to the statistical study of the old NSs in binaries in 
the aim at exploring the recycling process and the possibility that ultra fastly 
spinning NSs can form in these systems. The evolutionary scheme is based on 
§2,3.3, and §4.1. A similar approach is presented in §6 for the isolated NSs in 
the Galaxy. 

On the observational side, the NS having the shortest rotational period 
-Pmin = 1-558 ms ever detected is PSR 1937+21. Despite its apparent small- 
ness, P m in is not a critical period for NS rotation: as shown by Cook, Shapiro 

6 Teukolsky J8(J, the period P m i n is longer than the limiting period, P sn , below 
which the star becomes unstable to mass shedding at its equator, irrespective to 
the adopted equation of state (EoS) for the nuclear matter. Table [l] just shows the 
values of P sn for a set of EoS and the corresponding values of accreted baryonic 
mass. It shows how sensitive is P sn to the EoS, indicating that its determination 
is important for our understanding of nuclear processes in dense matter. 

The other important processes which intervene in accelerating a NS up to 
P s h are the evolution of the mass transfer rate from the low mass companion 
star and the evolution of the NS magnetic moment /x . Within the recycling 
scenario J8^] all schemes suggested for the origin and the evolution of /i allow 
the accreting NSs to reduce their [i and P to the values that are characteris- 
tic of MSPs. However, due to the large volume of the parameters in the mass 
transfer scenario and to the strong coupling between it and the field evolution 
of the NS in a low mass binary (LMB), only through a statistical approach it 
is possible to establish how efficient the recycling process is in spinning a NS to 
P < Pmin- This was recognized first by Possenti et al. f84|[[3l], who carried on 
statistical analyses of NSs in the millisecond and sub-millisecond range, using 
a Monte Carlo population synthesis code with ~ 3000 particles. It accounts for 
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Table 1. Disc accretion onto a 1.4 Mq neutron star 



EoS 


Fate 




AM B 


Mr* f,-„ 


P f 




-^rotating 
G,max 


nabs 

min 


F 


Collapse 


0.034 


0.172 


1.52 


0.72 


1.46 


1.67 


0.47 


A 


MassShed 


0.120 


0.428 


1.77 


0.60 


1.66 


1.95 


0.47 


E 


MassShed 


0.115 


0.414 


1.76 


0.66 


1.75 


2.05 


0.48 


ATT 
AU 


MassShed 


0.126 


0.446 


1.79 


0.70 


2.13 


2.55 


0.47 


D 


MassShed 


0.111 


0.405 


1.76 


0.73 


1.65 


1.95 


0.57 


FPS 


MassShed 


0.114 


0.416 


1.76 


0.75 


1.80 


2.12 


0.53 


UT 


MassShed 


0.119 


0.429 


1.78 


0.75 


1.84 


2.19 


0.54 


UU 


MassShed 


0.121 


0.436 


1.78 


0.78 


2.20 


2.61 


0.50 


C 


MassShed 


0.103 


0.389 


1.74 


0.89 


1.86 


2.17 


0.59 


N* 


MassShed 


0.130 


0.484 


1.84 


1.08 


2.64 


3.22 


0.68 


L 


MassShed 


0.116 


0.443 


1.80 


1.25 


2.70 


3.27 


0.76 


M 


MassShed 


0.091 


0.367 


1.74 


1.49 


1.80 


2.10 


0.81 



The second and the third columns contain the fate of the NS at the end point of re- 
cycling and the resulting ratio of kinetic energy over gravitational energy. The fourth 
and fifth columns report the total accreted baryonic mass and the final gravitational 
mass (in units of solar masses) . The sixth column lists the final attained rotational pe- 
riods. The seventh column collects the values of the maximum mass for a non-rotating 
spherical configuration, the eighth for a maximally rotating star and its correspond- 
ing minimum period for stability (ninth column). The table is from Cook, Shapiro & 
Teukolsky |o). 

the evolution during the early phases when the NS in the LMB behaves as if 
isolated (£ and V) and later when fed by wind accretion (A). The mass transfer 
during RLO (again A) is modelled considering a range of accretion rates which 
is close to the one observed. The population synthesis code follows also the last 
radio-pulsar phase, whose duration is chosen from a flat probability distribution 
in the logarithm of time (see Table ^ for a summary of the parameters used 
in the population synthesis code). The model incorporates the detailed physics 
of the evolution of a crustal magnetic field (as discussed in §3.3, using BCI and 
BCII to mimic expulsion or assimilation of the field in the NS core), and includes 
the relativistic corrections j8l[] necessary to describe the spin-up process. 

Evolution is followed also beyond the RLO A phase when accretion termi- 
nates. The increasing evidence that NSs in LMBs may suffer phases of transient 
accretion (perhaps due to thermal-viscous instabilities in an irradiation dom- 
inated disc ^5||8(|) is suggestive that mass transfer onto a NS may not stop 
suddenly: the star probably undergoes a progressive reduction of the mean ac- 
cretion rate, modulated by phases of higher and lower accretion. This in turn 
can start a cycle of V phases which in principle could vanish the effect of the 
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Table 2. Population syntheses parameters 



Physical quantity 



Distribution 



Values 



Units 



NS period at t£ LO (*) 

NS fj, at t£ LO (*) 
M in RLO phase (") 
Minimum accreted mass 
RLO accretion phase time (') 
MSP phase time 



Flat 
Gaussian 
Gaussian 
One- value 
Flat in Log 
Flat in Log 



1 -> 100 sec 

Log< no >=28.50 ; a=0.32 G cm 3 

Log< M >= 1.00 ; <r=0.50 Me 

0.01 M 

10 6 - T$Eh(*) year 

10 s -> 3 x 10 9 year 



(*) to = initial time of the Roche Lobe Overflow phase 

(jt) baryonic accretion rate during the Roche Lobe Overflow phase 

(f) a Maximum accreted Mass of 0.5Mq is permitted during the RLO phase 

(J) max duration of the RLO phase; explored values: 5xlO r yr - 10 s yr - 5xl0 8 yr 



previous spin-up. With the aim of exploring the effect of a decaying M on the 
population of fastly spinning objects, two possibilities have been investigated: a 
persistent accretion for a time trlo o r a persistent accretion for a shorter time 
followed by a transient phase mimicking the quenching of the mass transfer. The 
quenching of accretion has been modelled as a power law decay for M with in- 
dex P varying from 1 — > 10 (the last value is representative of an almost sudden 
switch off). 

Figure || collects the fractional distribution of the recycled model NSs using 
the set of parameters reported in Table || with TrIq = 5 x 10 8 yr. Guided by the 
values of P m in and /x m i n (the shortest rotational period observed in PSR 1937+21 
and the weakest magnetic moment observed in PSRJ2317+1439), the particles 
are divided in four groups. Those filling the first quadrant (P > P m i n and fi > 
A*min) behave as the known MSPs. Also the objects belonging to the second 
quadrant (P < P m j n and \i > fi m i n ) should shine as PSRs |8j. The effective 
observability of the objects in the third quadrant (P < P mm and /i < fi m i n ) 
as radio sources represents instead a challenge for the modern pulsar surveys. 
Most of them will be above the "death- line" j8g] , and might have a bolometric 
luminosity comparable to that of the known MSPs. Thereafter we shortly refer 
to su6-MSPs as to all objects having P < P m - m and \i above the "death-line". 
Objects in the fourth quadrant (P > P m - m and /i < fi min ) are probably radio 
quiet NSs, because they tend to be closer to the theoretical "death- line" , and 
they are in a period range which was already searched with good sensitivity by 
the radio surveys. 

The four crosses (Fig. ^) are for two representative EoSs (a very stiff EoS 
with a break-up period 1.4 ms and a mildly soft one, with P s h — 0.7 ms) and for 
two different boundary conditions (BCI & BCII) for the magnetic field at the 
crust-core interface. It appears that objects with periods P < P m i n are present 
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Fig. 6. Distributions of the synthesized NSs, derived normalizing the sample to the 
total number of model stars with P < 10 ms and arbitrary value of /x. The fj, — P plane 
is divided in four regions. As a guideline the upper left number in each cross gives the 
percentage of objects having P < P m in and /i > /x m in (the typical variance is about 
1%). 



in a statistically significant number. In effect, a tail of potential sub-MSPs always 
emerges for any reliable choice of the parameters listed in Table Q 

In the synthetic sample for the mild-soft EoS, the "barrier" at P sri ~ 0.7 
ms is clearly visible in Figure |^ (solid lines): the mild-soft EoS gives rise to 
period- distributions that increase rather steeply toward values smaller than 2 ms, 
irrespective to the adopted BCs. Instead, the boundary condition affects the 
distribution of /i: BCII produces a smaller number of objects with low field. 
Thus, sm6-MSPs could be a tool to test the physics of magnetic field decay. | 

Even the very stiff EoS permits periods P < P m i n , but the "barrier" of mass 
shedding (at P s h — 1.4 ms) is so close to P m - m = 1.558 ms that only few NSs 
reach these extreme rotational rates. Moreover the period distribution for the 
stiff EoS is much flatter than that for the soft-EoS, displaying a broad maximum 
at P ~ 3 ms. It has been recently claimed that X-ray sources in LMBs show 
rotational periods clustering in the interval 2^4 ms |sj| |x| . This effect could 
be explained introducing a fine tuned relation between \i and M(u oc M 1 / 2 
| fj9| ) . Alternatively, gravitational wave emission has been invoked |)l| |92) . Here 
we notice that such a clustering can be a natural statistical outcome of the 
recycling process if the EoS for the nuclear matter is stiff enough. 

Physical ingredients to describe the propeller induced spin-down of a NS at 
the end of the RLO phase are poorly known or difficult to assess (e.g. the exact 
law for the decrease of the mass transfer rate or the efficiency in the extraction 
of the angular momentum from the NS to the propelled matter). In our scheme, 



3 Note that a simple spin-down induced decay scenario has difficulty in explaining the 
observed MSPs and predicts no ultra fastly spinning objects, as shown in |84]. 
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Fig. 7. Calculated distribution of millisecond NSs as a function of the spin period P 
(fj, is let vary over the whole range). Solid line denotes the distribution in absence of 
propeller, whilst dashed area denotes the case of a strong propeller effect. The absolute 
number of objects is in arbitrary units. 



the largest effect occurs when phase V lasts half (~ 50%) of the RLO phase and 
the power-law index is T ~ 8. Figure reports the period distributions (dashed 
areas) when such a strong propeller phase is included. We note that a strong 
propeller can threaten the formation of NSs with P < P m i n and [i > /i m i n in the 
case of the very stiff EoS, whilst for the mild-soft EoS the distributions preserve 
a maximum just about P m m- 

This statistical analysis provides also information on the NS mass distribution 
as a function of P at the end of evolution. The initial NS gravitational mass, 
in the evolutionary code, was set at Mq = 1-35 M Q (according to the narrow 
Gaussian distribution with a = 0.04 M Q resulting from measuring the mass of 
the NSs in five relativistic NS-NS binary systems p3|]). Figure || clearly shows 
that the observed millisecond population should have undergone a mass load 
of ,$ O.1M during recycling. This is consistent with the few estimates of the 
masses of millisecond pulsars in low mass binaries |m| and raises the problem 
of explaining how the NS can get rid of the bulk of the mass (0.5 — 1.5 Mq) 
released by the companion during the RLO phase Jf2| . 

Figure || suggests that the mass function steepens toward high values only 
when P falls below ~ P m in, approaching M ~ 1.7 -j- 1.8 M a . That is a straight 
consequence of a results already pointed out by Burderi et al. jsjj : a large mass 
deposition (at least ^ 0.25M Q ) is required to spin a NS to ultra short periods, as 
illustrated in Figure |9[ The action of the propeller during the evolution has the 
effect of reducing the mass infall: the mass distribution is only slightly affected 
for the mild-soft EoS, while for the stiff EoS the difference is more pronounced. 
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Fig. 8. Average gravitational masses of the re-accelerated NSs as a function of their 
final spin period P. jj, varies over the whole range and the synthesized NSs are binned 
in 0.5 ms wide intervals. The initial mass of the static NSs is set equal to 1.35M in 
all the cases. Thick solid line denotes the mass-distribution in absence of propeller, 
whereas thin dashed line with a strong propeller included. 
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Fig. 9. Period versus Accreted Baryonic Mass for a magnetized NS using FPS-EoS 
and an initial gravitational mass of 1.40 Mq. The different pathways (a sample is 
represented by the dashed lines) define a strip, which narrows towards shorter periods. 
The strip is upper bounded by the evolutionary path for an unmagnetized NS (bold 
solid line). The bold long dashed line refers to the evolution along the spin- up line and 
is calculated assuming a tuned torque function (which maximizes the efficiency of the 
spin- up process). The thin solid line refers to a very slow decay of fi, implying larger 
mass depositions in order to attain very short spin periods. 
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The wandering among the various phases can account for the interesting 
phenomenology of the MSPs observed and opens the possibility that even more 
extreme objects like the sm&-MSPs exist in the Galaxy. 



6 The NS census in the Milky Way 

The isolated NSs follow a completely different evolution pattern. Generally, over 
the Hubble time, they cover large portion of the Milky Way and thus explore 
regions where the ISM is inhomogeneous looping among the various phases er- 



ratically. With a population synthesis model, Popov et aL[101] traced their evo- 
lution in the ISM, and in the Milky Way potential. The model NSs, born in 
the galactic plane at a rate proportional to the square of the ISM density, have 
initially short spin periods, magnetic fields clustered around 10 30 G cm 3 , and 
spatial (kick) velocities that can be drawn for a Maxwellian distribution (with 
mean velocity modulus (V) treated as a parameter). 
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Fig. 10. Fractions of NSs (in percents) in phases E^VA and Q versus the mean kick 
velocity (V), for a constant magnetic moment /i = 0.5 x 10 30 Gcm 3 (open circles) and 
(i — 10 30 Gcm 3 (filled circles); typical statistical uncertainty for £ and A ~ 1-2%. 



The collective properties are illustrated in Figure |l0| for two initial values of 
/io = 0.5 — 1 x 10 30 G cm 3 : For a non evolving field most of the isolated NSs 
spend their lives as Ejectors and there is no possibility to observe them after 
the radio pulsar phase. Propellers are shortlived 0, and Georotators are rare. 
A tiny fraction (a few percents) on the NSs are in the Accretor stage if (V) is 
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above 200 km s 1 . Only in the unrealistic case of a low mean velocity, the bulk 
of the population would be in stage A. 

As illustrated in Figure O, in phase A, the mean value of M is a few 10 9 g 
s" 1 and, among the Accretors, the typical velocity clusters around 50 km s , 
and the luminosity L ~ (GM / R)M around 10 29 erg s _1 . Accreting isolated NSs 
are "visible", but they would be extremely dim objects emitting predominantly 
in the soft X-rays, with a polar cap black body effective temperature of 0.6 keV. 
What can we learn when comparing the theoretical census with the observations? 
We can discover if field decays in isolated objects, if the population is devoid 
of slow objects, and if the stars effectively spin down due to the unavoidable 
interaction with the tenu ous ISM. The search of accreting isolated NSs is thus 
compelling |(j§ Eel |og Q . 




log S 
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Fig. 11. Upper left panel: the log N-log S distribution for Accretors within 5 kpc from 
the Sun. Dashed (solid) curve refers to emission from polar cap (entire NS surface) in 
the range 0.5-2.4 keV: straight lines with slopes -1 and -3/2 are included for comparison. 
Fro top right to bottom right: the velocity V, effective polar cap temperature and 
accretion rate M distributions for Accretors [106 



6.1 Accreting isolated neutron stars in the Rosat Sky? 

Despite intensive observational campaigns, no irrefutable identification of an 
isolated accreting NS has been presented so far. Six soft sources have been found 
in Rosat field j78|, identified as isolated NSs from the optical and X-ray data. 
Observations, however, do not permit to unveil the origin of their emission, yet. 
These sources could be powered either by accretion or by the release of internal 
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energy in relatively young (as 10 6 yr) cooling NSs [jj. Although relatively bright 
(up to « 1 ctss -1 ), the proximity of the sources (inferred from the column 
density in the X-ray spectra) makes their luminosity (L « 10 31 ergs -1 ) near to 
that expected from either a close- by cooling NS or from an accreting NS, among 
the most luminous. Their X-ray spectrum is soft and thermal, again as predicted 
for both Accretors and Coolers. 

Can a comparison with theoretical expectations help in discriminating among 
the two hypothesis? First, the paucity of very soft X-ray sources in the Rosat 



fields (in comparison with earlier expectation [103] [ 104 1) is indicating that Ac- 
cretors are rare objects. If these six sources are indeed accreting, this implies 
that the average velocity of the isolated NSs population at birth has to exceed 
~ 20 kms -1 (a figure consistent with that derived from radio pulsars statistics 
[ [105| ). In addition, since observable accretion-powered isolated NSs are (intrinsi- 
cally) slow objects, these results exclude the possibility that the present velocity 
distribution of NSs is richer in low-velocity objects with respect to a Maxwellian. 
Thus, despite the fact that in our Galaxy there are many old NSs (about 10 9 ), 
the young (~ 10 6 yr) cooling NSs seem to outnumber, at the Rosat counts, those 
in phase A. This is what emerges also from the calculation of the log N-log S dis- 
tribution both of cooling and accreting stars (the last carried on with " census" ; 
[ |106| ) . At the bright counts, the local population of cooling NSs would dominate 
over the log N-log S of the dimmer and warmer (less absorbed) Accretors. There 
is the hope that Chandra and Newton will detect them, at the flux limit of 10 -16 
erg cm~ 2 s _1 106( . In support of the cooling hypothesis there is also the recent 
measurement [107 1 [108 of the velocity (200 km s _1 ) of RXJ1856 (a member of 
this class) indicating that this source is, most likely, not powered by accretion 
having a high velocity. Interestingly, the cooling hypothesis (when explored in 
the log N-log S plane) implies that the NS birth rate in the Solar vicinity, over 
the last 10 6 yr, is higher than that inferred from radiopulsar observations [106 



[ |109| . This might be a crucial point as it may indicate that most of the young 
NSs may not share the observational properties of the canonical PSRs, hinting 
for the presence of a background of "anomalous" young cooling NSs. 

What about field decay in the accretion hypothesis? Decay of the field to 



extreme low values would imply a large number of Accretors [ 10 1 1 , not observed, 
thus excluding this possibility. This indicates that, in crustal models Q |}(J, 
electric currents need to be located deep into the crust and that the impurity 
parameter needs to be not exceedingly large |j4| , in spin-down induced models. 
What is still uncertain is whether the paucity is due mainly to a velocity effect 
than to a "fine-tuned" decay, a problem not solved yet, statistically. 

While our theoretical expectation hints in favor of the cooling hypothesis 
there remain yet a puzzle: the long period of one of these sources, RXJ0720 
floj . If powered by accretion, RXJ0720 would have a "weak" field NS {p, > 10 26 
G cm 3 ) and this would be a rather direct prove of some field decay, in isolated 



NSs|q] |lll| . Can a young cooling object have such a slow rotation? Do we have 



see 178| for an updated review, the description of the sources and a complete reference 
list. 
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to change or view that NSs come to birth with ultra short periods? A new 
challenging hypothesis has been put forward [112 |113| that RXJ0720 is 
just the descendant of a highly magnetized NS (a Magnetar) that have suffered 
a severe spin down accompanied by a nonlinear decay of the field whose energy is 
powering the X-ray luminosity 114]. This issue remains one of the new problems 
of the NS physics, making the debate on the nature of these sources an even more 
exciting problem. 



7 Conclusions 



In this review we have traced the evolution of old NSs transiting through the 
Ejector, Propeller and Accretor phases. The NSs of our review are far from 
being "canonical" MSPs in light binaries, or "canonical" PSRs in the field. The 
ultra fastly spinning NSs, that we have recycled in binaries, are rather extreme 
relativistic heavy NSs: If discovered, they will unable us to probe the stellar 
interior in an unprecedented way, and to constrain the physics driving magnetic 
field decay, in interacting systems. As regard to the field NSs, the discovery of the 
six Rosat sources has just opened the possibility of unveiling "unconventional" 
NSs, evolving in isolation. Whether they are accreting or cooling objects is still 
a mystery and even more mysterious and fascinating is their possible link with 
" Magnetars" . The study of these " unusual" NSs can open new frontiers in this 
already active field of research. 
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